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PREFACE 


This  pamphlet,  one  cf  the  series  which  will  comprise  the  Ordnance  Engineering 
Design  Handbook,  deals  with  the  subject  cf  Exterior  Ballistics  ,  covering  Tra¬ 
jectories,  Differential  Effects  and  Data  for  Projectiles.  The  pamphlet  pertains 
to  that  phase  of  projectile  life  starting  at  the  gun  muzzle  (or,  for  a  rocket,  at 
burnout)  and  ending  at  the  point  of  burst  or  impact.  The  subject  is  covered  in  a 
manner  intended  primarily  to  assist  the  ammunition  designer  in  recognizing  and 
dealing  with  the  design  parameters  of  projectile  weight,  shape,  muzzle  velocity,  yaw, 
drag  and  stability,  in  designing  for  optimum  results  in  range,  time  cf  flight  and 
accuracy.  It  is  hoped  that  it  will  also  prove  helpful  to  other  personnel  cf  the  Ord¬ 
nance  Corps  and  of  contractors  to  the  Ordnance  Corps.  Because  of  its  complexity, 
it  is  possible  to  cover  the  subject  only  briefly  in  this  volume,  but  an  attempt  has  been 
made  to  incorporate  references  which  will  permit  the  designer  or  the  student  to 
explore  in  more  detail  any  phase  of  the  subject. 

The  ballistics  cf  bombs  and  of  intennediate  and  long  range  missiles  are  not 
covered  in  this  pamphlet.  Neither  are  the  problems  associated  with  Interior 
Ballistics  or  of  Terminal  Ballistics.  These  subjects  will  be  covered  by  other 
pamphlets. 
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SYMBOLS 


A — Azimuth,  measured  from  the  direction  of  mo¬ 
tion  of  the  airplane  (in  aircraft  fire) 

A — Altitude  function  (Eq.  63) 

A — Axial  moment  of  inertia 
B — Drag  coefficient  in  lb/in2-ft  (Eq.  19) 

C — Ballistic  coefficient  (Eq.  22) 

D — Drag 

D— A  ngular  deflection 
F  —  Temperature  in  degrees  Fahrenheit 
G — Drag  function  (Eq.  20) 

H — Ratio  cf  air  density  at  altitude  y  to  that  at 
the  surface  (Eq.  48) 

H — Height  of  target 
I — Inclination  function  (Eq.  62) 

K — Coefficient  (See  Eq.  16  and  34) 

L — Crosswind  force  (lift) 

M — Mach  number  (Eq.  18) 

M — Pseudo-Mach  number  (Eq.  65) 

M  —  Overturning  moment  about  the  center  of 
gravity 

M — log10e=0.43429  (approx.) 

N — Spin 

P  —  Performance  parameter  (Eq.  27) 

P — Siacci  range  :  the  distance  measured  along  the 
line  of  departure  to  a  point  directly  above 
the  projectile  (Eq.  54) 

Q — Drop  of  the  projectile  (Eq.  55) 

R — Reynolds  number  (Eq.  17) 

R — Range 

S — Space  function  (Eq.  60) 

T  —  Total  time  cf  flight 
T  —  Time  function  (Eq.  61) 

U — Upper  limit  d  integration  of  Siacci  functions 
V — Velocity 
W — Wind  speed 
X — Total  range  cf  a  trajectory 
Y  —  Altitude  d  an  airplane  above  sea  level 

iv 


2 — Linear  deflection 
2 — Zenith  angle 
i — Velocity  of  sound 

Diameter  d  the  largest  cross  section  of  the 
projectile  (in  Givre  drag  law,  Eq.  23) 
o' — Damping  coefficient  due  to  yawing  moment 
and  crosswind  force  (Eq.  51) 
c" — Damping  coefficient  due  to  drag  (Eq.  33) 
d — Caliber,  or  maximum  diameter  of  projectile 

f — Drag  coefficient  (in  the  Givre  drag  law,  F.q. 
23) 

g — Gravitational  acceleration 
h — Logarithmic  rate  of  decrease  of  air  density 
with  altitude:  0.000,031,58  per  ft 
h — Yawing  moment  damping  factor 
i — Form  factor  (Eq.  21) 

k — Logarithmic  rate  of  decrease  of  sonic  velocity 
with  altitude  0.000,003,01  per  ft 
k — Ratio  cf  American  to  French  drag  function 
(Eq.  24) 

m — Mass  of  the  projectile 
m — Slope  of  the  tangent  to  the  trajectory 
n-4  .ogarithmic  rate  if  change  of  muzzle  velocity 
with  projectile  weight 
p — Space  function  (Eq.  28) 
p — Projectile  weight 
qi — Inclination  function  (Eq.  30) 
q — Altitude  function  (Eq.  31) 
s — Stability  factor 

t — Time  of  flight  (sometimes  the  symbol  tf  is  used 
to  distinguish  the  time  d  flight  from  the 
Siacci  tinie  function) 
t — Time  function  (Eq.  29) 
u — Velocity  of  the  projectile  relative  to  the  air 
u — Component  along  the  line  d'  departure  of  the 
velocity  relative  to  the  air:  the  argument  of 
the  Siacci  functions 

v — Velocity  of  the  projectile  relative  to  an  inertial 
system 


w — True  air  speed  cf  the  airplane  in  aircraft  fire 
x  —  Horizontal  distance  from  the  origin  cf  the 
trajectory 

y — Vertical  distance  from  the  origin  of  the  tra¬ 
jectory 

z — Drift 

A  (capital  delta)  or  A' — Increment  cf  (fbUctrad. 
by  a  variable) 

A — Density  cf  air  (in  Gavre  drag  law,  By.  23) 
y  (gamma) —  Factor  inversely  proportional  to  the 
ballistic  coefficient  (Eq.  66) 

3  (delta) — Yaw 

4  (epsilon) — Anglecf  site  (Eq.  72) 

iS  (theta)  —  Angle  of  inclination  cf  the  trajectory 
%  (kappa) — Crosswind  force  damping  factor 
|i  (mu)  —  Viscosity  cf  the  air 
p  (rho)  —  Density  cf  the  air 
t  (tau)  —  Temperature  cf  the  air 
f  (phi)  —  Angle  of  departure 
q  (omega)  —  Angle  cf  fall 
Superscripts 

One  dot  (example  :  £)• — Derivative  with  respect 
to  time. 

Two  dots  (example  :  y) — Second  derivative  vriiii 

respect  to  time. 

Subscripts 

0— At  t  =  0,  y  =  0,  or  3  =  0, 


1  — Pertaining  to,  or  relative  to.  Projectile  Type  1 
2 — Pertaining  to,  or  relative  to,  Projectile  T^pe  2 
8 — Pertaining  to,  or  relative  to.  Projectile  Tgpe  8 
.1  — First  revision 
.2 — Second  revision 

C — Due  to  1%  increase  in  ballistic  coefficient 

D — Drag 

H  — Due  to  height  cf  target 
R — Range  component  cf 
e — Effective 

p — Due  to  projectile  weight 
s — Summital  :  at  the  summit  cf  a  trajectory,  where 

y  —  o 

s — Standard  :  pertaining  to  standard  conditions 
t  —  Pertaining  to,  or  relative  to,  some  Typical  Pro¬ 
jectile 

v — Due  to  1  unit  increase  in  muzzle  velocity 
w — Due  to  wind 
z — Cross  component  cf 

5— Yaw 

t — Due  to  increase  in  the  temperature  cf  the 
atmosphere 

? — Due  to  an  increase  in  angle  cf  departure 

1 Terminal:  where  y  =  0  on  the  descending 
branch  cf  the  trajectory 


TRAJECTORIES,  DIFFERENTIAL  EFFECTS,  AND  DATA 

FOR  PROJECTILES 


1.  Introduction  to  Extejuor  Ballistics 

Ballistics  is  the  science  cf  the  motion  and 
effects  of  projectiles.  Exterior  Ballistics  deals 
with  that  part  of  the  motion  between  the  pro¬ 
jector  and  the  point  cf  burst  or  impact. 

The  projector  may  be  any  sort  cf  weapon:  a 
sling-shot,  a  bow,  a  mortar,  a  rifle,  a  rocket 
launcher,  an  airplane,  etc.  A  projectile  may  be 
any  body  that  is  projected  by  such  a  weapon. 
However,  in  this  pamphlet,  the  term  “projectile” 
will  be  restricted’to  those  missiles  that  are  in  free 
flight  through  the  air  after  being  propelled  by  a 
gun  or  a  rocket  motor.  This  term  does  not  in¬ 
clude  bombs,  which  are  released  from  airplanes; 
they  will  be  considered  in  another  pamphlet. 

The  motion  of  a  missile  in  flight  is  really 
quite  complicated.  Even  a  fin-stabilized  pro¬ 
jectile  or  bomb  that  is  not  spinning  about  its 
longitudinal  axis  is  swinging  about  a  transverse 
axis,  and  this  affects  its  motion.  The  nose  of  a 
spinning  shell  precesses  about  the  tangent  to  its 
trajectory,  and  its  attitude  gives  rise  to  forces  that 
not  only  affect  its  motion  in  a  vertical  plane  but 
also  move  it  out  cf  this  plane  To  simplify  the 
calculations,  however,  it  is  usually  assumed  that 
a  missile  acts  like  a  particle  subject  to  two  forces : 
(1)  gravity,  which  is  directed  vertically  down¬ 
ward,  and  (2)  drag,  which  is  directed  opposite  to 
the  direction  of  motion.  Gteaaty  is  assumed  con¬ 
stant  :until  1956, the  gravitational  acceleration  was 
usually  taken  to  be  32.152  ft/sec*  (S.SOm/sec4)  ; 
and  was  used  for  all  tables  referred  to  in  this  text; 
now,  the  standard  value  is  32.174  ft/sec*  (9.80665 
m/sec5).1  The  drag  depends  on  the  size  of  the 
missile,  its  velocity,  the  density  <f  air,  and  other 
factors.  Under  these  conditions,  the  motion  fol¬ 
lows  a  smooth,  planar  path,  which  is  called  an 
“exterior  ballistic  trajectory,”  or  “trajectory”  for 
short 

The  initial  conditions  for  the  trajectory  cf  a 
particular  missile  are  its  initial  velocity  and  angle 
cf  departure.*  For  a  projectile,  the  initial  velocity 

*  The  term  attgk  of  departure  and  the  term  elevation 
as  used  in  this  pamphlet  represent  the  angle  from  the 
korisontal  to  the  line  cf  departure  and  therefore  would  be 
more  accurately  and  completely  termed  quadrant  angle  of 
departure. 


is  the  muzzle  velocity  and  angle  of  departure  is 
approximately  the  elevation  cf  the  weapon 
(fictitious  values  of  these  quantities  may  be  used 
in  the  case  cf  a  rocket  that  bums  outside  the 
launcher).  A  bomb  is  released  with  the  speed  of 
the  airplane,  generally  in  a  horizontal  direction. 
If  a  gun  is  mounted  in  an  airplane  or  other 
moving  vehicle,  the  initial  conditions  are  affected 
by  the  motion  cf  the  vehicle.  In  any  case,  it  is 
customary  to  compute  trajectories  for  a  few  se¬ 
lected  values  cf  initial  velocity  and  elevation. 

Certain  standard  conditions  are  specified  for  a 
normal  trajectory  :not  only  the  initial  velocity  and 
angle  cf  departure,  but  also  the  weight  of  the 
projectile  and  a  factor  that  depends  on  its  shape: 
the  density,  temperature,  and  lack  of  motion  cf 
the  air,  and  the  height  of  the  target  relative  to  the 
projector.  In  order  to  detennine  the  ‘effects  cf 
variations  from  these  conditions,  it  may  be  con¬ 
venient  in  some  cases  to  compute  non-standard 
trajectories;  however,  they  can  usually  be  esti¬ 
mated  in  other  ways,  as  will  be  explained  later. 

2  Vacuum  Trajectories 

a.  Equations  of  Motion 

A  rough  idea  of  an  exterior  ballistic  trajectory 
may  be  gamed  by  neglecting  the  air  resistance, 
considering  only  the  effects  of  gravity.  Then  if 

x  is  the  horizontal  distance  from  the  origin, 
y  is  the  vertical  distance  from  the  origin, 
t  is  the  time  cf  flight, 
g  is  the  gravitational  acceleration, 

and  dots  denote  derivatives  w'ith  respect  to  time, 
the  differential  equations  of  motion  are 

x  =  0,  y  —  —  8-  ( 1 ) 

The  initial  ''onditions,  denoted  by  the  subscript  0, 

are 

x*  =  v0  cos  ■$<,,  y„  =  v0  sin  =  0,  y0  =  0, 

where  v  is  the  velocity, 

■S  the  angle  cf  inclination. 

By  integrating  equations  (1),  we  find  that,  at 
any  time!, 

x  =  v0  cos-#*,  y  —  v„  sir;  —  gt,  (2) 

and  by  integrating  (2) , 
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x  =  (v0  cos  d0)t,  y  =  (v0  sind0)t  — gt s/2.  (3) 
b.  Characteristics  of  Trajectories 

From  Equations  (3),  we  find  that  the  time  at 
any  horizontal  range  is 

t  =  X/Vo  cos  d0  (4) 

and  hence  the  ordinate  is 
y  =  x  tan  -So  —  gx2/2v02  cos2  •$„.  (5) 

This  is  a  parabola  with  a  vertical  axis. 

The  total  range  X  is  found  by  letting  y  =0  in 
(5)  and  solving  for  x: 

X  =  (2v„2/g)  sin£0cos£0  =  (v02/g)  sin2*0.  (6) 

This  shows  that  the  range  is  a  maximum  when 
a,  =  45". 

The  total  time  of  flight  T  is  found  by  substi¬ 
tuting  (6)  for  x  in  (4)  : 

T  =  (2v0/g)  sin$o-  (7) 

The  components  of  the  terminal  velocity  are 

found  by  substituting  (7)  for  t  in  (2)  : 

X  =  v0cos£O)  Y  = — v0  sind0.'  (8) 

Hence  the  terminal  velocity  is 
v<„=v0,  (9) 

and  the  angle  of  fall  (considered  positive)  is 
<a  =  40.  (10) 

By  differentiating  equation  (5)  with  respect  to 
x,  we  find  that  the  slope  of  the  trajectory  at  any 
point  is 

=  tan  —  gx/v02  cos2  d0.  (11) 

At  the  summit,  the  slope  is  0;  hence  the  horizontal 
range  to  the  summit  is 

x,  =  (v„*/g)  sin  £0  cos  $0  =  X/2.  (12) 

By  substituting  this  for  x  in  (5),  we  find  the  max¬ 
imum  ordinate 

y.  =  (v02/2g)  sin2  d0  =  (1  —  cos230)  v0!/4g 

—  (g/8)  T3.  (13) 

By  substituting  (12)  for  x  in  (4),  we  find  the 

time  to  the  summit 

t.  =  (v0/g)  sin  .=  T/2.  (14) 

3.  Air  Resistance 

a.  Forces  Acting  on  a  Missile  Moving  in  AV 
A  missile  moving  in  air  encounters  other  forces 
besides  gravity.  The  principal  components  of  the 
force  produced  by  motion  through  air  are  drag, 
crosswind  force,  and,  for  a  spinning  projectile, 
Magnus  force 

( 1 )  Drug  is  the  component  in  the  direction  op¬ 


posite  to  that  c£  the  monon  of  the  center  of 
gravity.  This  is  the  only  component  acting  on  a 
non-spinning  projectile  that  is  trailing  perfectly. 

(2)  Crosswind  Force  in  the  component  per¬ 
pendicular  to  the  direction  of  motion  of  the  center 
of  gravity  in  the  plane  cf  yaw.  The  yaw  is  the 
angle  between  the  direction  of  motion  cf  the 
center  cf  gravity  and  the  longitudinal  axis  cf  the 
missile.  This  component  is  the  cause  cf  the  drift 
cf  a  spinning  shell. 

(3)  Magnus  Force  acts  in  a  direction  per¬ 
pendicular-  to  the  plane  of  yaw.  This  is  the  force 
that  makes  a  spinning  golf  ball  or  baseball  swerve. 

b.  Simplified  Assumptions 

The  motion  caused  by  these  forces  is  explained 
in  detail  in  the  Ordnance  Corps  Pamphlet,  “De¬ 
sign  for  Control  of  Flight  Characteristics.”2  In 
general,  the  motion  of  the  center  of  gravity  is 
affected  by  the  motion  about  the  center  of  gravity 
and  occurs  in  three  dimensions.  To  simplify  the 
calculations,  it  is  usually  assumed  that  the  only 
forces  acting  on  the  missile  in  free  flight  are 
drag  and  gravity.  The  crosswind  force  is  implied 
in  the  drift,  but  this  is  determined  experimentally 
as  a  variation  from  the  plane  trajectory.  This 
approximation  is  satisfactory  except  at  extremely 
high  elevations,  where  the  summital  yaw  is  very 
large. 

c.  Drag  Coefficient 

The  drag  encountered  by  a  projectile  moving  in 
air  is  defined  by  the  equation  cf  motion 
dv  . 

m  =  — D  —  mg  sin  *,  (15) 

where 

m— is  the  mass  cf  the  projectile, 
v — tiie  velocity  cf  the  projectile  relative  to  an 
inertial  system, 
t — tiie  time, 

D— tiie  drag, 

g — tiie  gravitational  acceleration, 

•S — the  angle  of  inclination. 

The  drug  coefficient  is  defined  by  the  relation 

KD  =  D/pd*u*,  (16) 

where 

KD — is  the  drag  coefficient, 
p— tiis  air  density, 
d — the  caliber  or  maximum  diameter, 
u — tiie  velocity  cf  the  projectile  relative  to 
the  air. 

The  drag  coefficient  is  a  function  cf  the  Reynolds 
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number  R  and  the  Mach  number  M.  The 
Reynolds  number  is  defined  as 

R=udp/|i  (17) 

where  jj.  is  tlie  viscosity  cf  tlie  air,  and  the  Efadi 
number  as 

M  =  u/a,  (18) 

where  a  is  the  velocity  of  sound  in  the  surround¬ 
ing  medium.*  The  viscosity  is  related  to  the 
friction  of  the  air;  the  velocity  cf  sound,  to  its 
elasticity.  At  extremely  low  velocities,  the 
Reynolds  number  causes  greater  variation,  but 
the  effect  of  the  Mach  number  is  larger  when  it 
exceeds  O.S.  The  angle  between  the  direction  of 
motion  of  the  center  of  gravity  and  the  axis  cf 
symmetry  cf  the  shell,  called  yaw,  also  affects  the 
drag  coefficient,  but  this  is  usually  stripped  out  cf 
the  resistance  firing  data. 

Figure  2  (p.  22)  is  a  graph  cf  drag  coefficient 
vs  Mach  number  for  the  four  typical  projectiles 
whose  shapes  are  shown  in  Figure  1  (p.  21). 
These  projectiles  are  described  and  the  basic 
firings  explained  in  Section  4.  It  should  be  noted 
that  subscripts  are  used  to  distinguish  the  typical 
projectiles  to  which  the  drag  coefficients  pertain; 
a  second  subscript  is  sometimes  added  after  a 
decimal  point  to  denote  a  revision.  It  is  evident 
that  all  the  drag  coefficients  increase  rapidly  in 
the  vicinity  of  Mach  number  1.  The  main  differ¬ 
ence  between  curves  is  in  the  ratio  of  the  maxi¬ 
mum,  which  occurs  between  Mach  numbers  1.0 
and  1.5,  to  the  minimum  subsonic  value.  Although 
the  drag  coefficient  has  been  determined  and  tab¬ 
ulated  for  several  other  projectiles,  these  four 
represent  the  principal  shapes  for  artillery  pro- 
j  ectiles. 

Kd  is  a  nondimensional  coefficient,  and  a  con¬ 
sistent  set  of  units  should  be  used  in  Equation 
(16).  However,  it  is  customary  to  express  the 
density  as  a  ratio,  the  caliber  in  inches,  and  the 
velocity  in  feet  per  second;  then,  since  the  drag 
is  in  poundals,  the  drag  coefficient  is  expressed  in 
pounds  per  square  inch  per  foot  and  is  denoted  by 
B.  Until  1956,  the  standard  air  density  was 
taken  as 

p0  =  0.07513  Ib/ft*  =  5.217  X  1(>-4  Ib/inMt,* 

The  practical  drag  coefficient, 

B  =  5.217  X  10-4  Kd,  (19) 

for  Projectile  'types  1,  2  and  8  is  tabulated  as  a 
function  cf  velocity  for  the  standard  velocity  <f 
sound, 

ao_=  1120.27  fps,+ 

*  In  1956,  the  Department  of  Defense  adopted  the  stand¬ 
ard  atmosphere  of  the  International  Civil  Aviation  Organ- 
hz±ki l'  In  this  atmosphere,  at  the  surface,  the  density  is 


The  nondimensional  drag  coefficient  Kd  for  the 
90-mm  HEAT  Shell  T108  is  tabulated  as  a  func¬ 
tion  cf  Mach  number.  (See  References  3,  4,  5  and 
6,  which  may  be  obtained  from  the  Ballistic  Re¬ 
search  Laboratories,  Aberdeen  Proving  Ground, 
Maryland). 

d.  Drag  Function 

The  drag  function  G  is  computed  by  the  formula 
G  =  Bu  (20) 

with  u  in  feet  per  second.  For  convenience  in 
computing  trajectories,  this  has  been  tabulated 
as  a  function  of  us/100  with  u  in  meters  per 
second  and  either  feet  or  yards  per  second,  for 
Projectile  Types  1,  2  and  8.  (See  References  7, 

8  and  9.) 

e.  Form  Factor 

It  is  not  practical  to  determine  the  drag  co¬ 
efficient  as  a  function  of  hjach  number  for  all 
projectiles.  Instead,  the  drag  coefficient  cf  a  pro¬ 
jectile  whose  shape  differs  only  slightly  from  one 
cf  the  typical  projectiles  is  assumed  to  be  pro¬ 
portional  to  the  typical  drag  coefficient.  If  KDj 
is  the  drag  coefficient  cf  a  Typical  Projectile,  the 
form  factor  of  a  projectile  whose  drag  coefficient 
is  Kd  is  the  ratio 

it  =  KD/KDt.  (21) 

The  form  factor  may  be  determined  from  the 
time  cf  flight  at  short  ranges,  or  from  the  range 
to  impact  or  the  position  cf  burst  at  a  few  eleva¬ 
tions.  If  it  is  desired  to  estimate  the  form  factor  cf 
a  projectile  without  firing  it,  the  Ordnance  Corps 
Pamphlet,  “Design  for  Control  of  Flight  Char¬ 
acteristics"2  should  be  consulted.  It  explains  the 
effects  on  drag  of  variations  in  length  and  radius 
cf  head,  meplat  diameter,  base  area,  length  of 
shell,  and  yaw. 

Generally,  the  form  factor  of  a  projectile  with 
a  square  base  and  an  ogive  less  than  1.75  calibers 
high  should  lie  referred  to  Projectile  Itype  1 ;  that 
cf  a  projectile  with  a  square  base  and  an  ogive 
more  than  1  75  calibers  high,  to  Projectile  Type 
8;  that  of  a  projectile  with  a  boattail,  to  Projectile 
■type  2;  and  that  of  a  finned  projectile,  to  the 
T108  Shell.  However,  there  are  exceptions  to 
these  rules  :e.g„  Kd,  is  used  for  fin-stabilized  mor¬ 
tar  shells  at  low  velocities,  and  Kr^  fits  the  calcu¬ 
lated  drag  coefficientfor  some  longer  fin-stabilized 
shells.  The  drag  coefficients  cf  bombs  are  treated 
elsewhere.10 

0.076,475  lb/ft*  and  the  velocity  of  sound  is  111689  fps; 
also  the  viscosity  is  1.205  x  10— *  lb/ft-sec,  but  tlie  effect 
of  Reynolds  number  is  not  considered  in  computing  tra¬ 
jectories. 
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/.  Ballistic  Coefficient 

The  ballistic  coefficient  relative  to  the  tabulated 
drag  coefficient  is  defined  by  the  formula 

Ct  =  m/itd2,  (22) 

where  m  is  the  mass  cf  the  projectile  in  pounds, 
d  is  the  caliber  (or  maximum  diameter) 
in  inches. 
g.  Air  Density 

The  air  density  is  a  function  of  altitude.  The 
exponential  function  adopted  by  the  Ordnance 
Corps*  is  tabulated  in  References  7  and  11.  The 
ratio  of  the  density  at  an  altitude  y  to  that  at  the 
surface  is  H(y). 

4,  Description  of  Typical  Projectiles 

This  pamphlet  contains  data  that  pertain  spe¬ 
cifically  to  four  types  of  projectiles:  Types  1,  2 
and  8  and  the  90-mm  HEAT  Shell  T108.  They 
can  be  applied  to  other  shapes  by  using  a  form 

factor,  as  explained  in  paragraph  3e. 

a.  Projectile  Type  1 

Most  of  the  service  projectiles  used  by  the 
European  artillery  during  the  last  third  cf  the 
nineteenth  century  had  square  bases  and  ogival 
heads  of  various  apical  angles.  In  Russia,  during 
the  years  1868  and  1869,  General  Magenski  con¬ 
ducted  some  experiments  at  St.  Petersburg  to  de¬ 
termine  the  law  of  air  resistance  of  projectiles  with 
a  semi-apical  angle  of  48°22'  (1.49-caliber  radius 
bead)  at  velocities  from  172  to  409  m/s.  In  Eng¬ 
land,  from  1866  to  1880,  Bashforth  conducted  sim¬ 
ilar  experiments  with  projectiles  cf  the  same  form 
(it  may  be  that  the  tables  were  reduced  to  this 
form)  at  velocities  from  148  to  825  m/s.  In  Hol¬ 
land,  Colonel  Hojel  conducted  a  few  such  experi¬ 
ments  during  1883;  the  exact  form  is  uncertain, 
but  the  ogive  was  evidently  somewhat  longer  than 
those  cf  the  Russian  experiments.  In  Germany, 
the  Krupp  factory  conducted  a  large  number  cf 
air  resistance  experiments  at  Meppen  Proving 
Ground,  starting  in  1881  (the  principal  results 
were  reported  by  Siacci  in.  the  "Rivista"  for 
March  18%)  ;  the  velocities  varied  from  368  to 
910  m/s,  and  three  different  types  of  projectiles 

*The  tabulated  trajectory  data  based  on  the  old  stand¬ 
ard  atmosphere  give  results  that  are  nearly  the  same  as 
those  based  on  the  new  atmosphere  providing  a  suitable 
change  is  made  in  the  ballistic  coefficient.  For  example,  a 
280-mm  shell,  fired  with  a  muzzle  velocity  of  1795  fps,  at 
an  elevation  of  45°,  attains  a  range  of  15,790  yds.;  the 
ballistic  coefficient  is  3.62  with  ICAO  atmosphere  and 
3.50  with  the  old  Ordnance  Corps  atmosphere.  This  func¬ 
tion  was  used  in  computing  trajectories  before  1956,  when 
the  Department  of  Defense  adopted  the  standard  atmos¬ 
phere  of  the  International  Civil  Aviation  Organization. 


were  -  Ve:  cr''  same  as  the  Russian  (or 

the  results  reduced  to  the  same),  but  the  others 
probably  had  longer  ogives.  In  France,  a  com¬ 
mission  cf  the  naval  artillery  conducted  many  such 
experiments  at  Gavre  Proving  Ground  from  1873 
to  1898  at  velocities  from  118  to  1163  m/s;  they 
fired  ogival  projectiles  with  three  semi-apical 
angles:  45°55'  (1.64  CRH),  41°40/  (1.98  CRH), 
and  31°45'  (3.34  CRH).  For  further  details,  see 
References  12  to  16. 

In  analyzing  the  results  cf  the  experiments,1* 
the  Givre  Commission  cf  Experiments  assumed 
that  the  drag  could  be  expressed  in  the  form 

D  =  a2V2  f(V),  (23) 

where 

D — is  the  drag  (kg), 

A — the  density  of  the  air  (kg/m*), 

g — the  gravitational  acceleration  (m/see1), 

a — die  diameter  of  the  largest  cross  section  of 
the  projectile  (m), 

V — die  velocity  (m/s). 

The  drag  was  treated  as  a  function  cf  velocity 
without  taking  account  of  the  velocity  cf  sound, 
although  the  dependence  cf  drag  on  the  velocity  cf 
soimd  was  strongly  suspected:  the  coefficient  f 
was  plotted  as  a  function  cf  V  alone.  After 
studying  the  results  obtained  with  projectiles  cf 
various  ogival  angles,  the  Commission  aenfirmad 
the  law  deduced  by  Helie  in  1888  that  f(V)  is 
proportional  to  the  sine  of  the  ogival  angle,  at 
least  for  values  from  40°  to  90°,  within  the  experi¬ 
mental  errors. 

In  1917,  Chief  Engineer  Gamier  cf  the  Givre 
Commission  published  a  table1*  cf  log  B  and  die 
corresponding  logarithmic  derivative,  with  ve¬ 
locity  as  argument:  this  B  is  equivalent  to  f(V) 
for  an  ogival  projectile  with  a  2-caliber  radius 
head.  This  ideal  shape  is  what  we  call  Projectile 
Tjjb  1  (see  Figure  1).  The  Ordnance 
ment  of  the  United  States  Army  then  prepared  a 
table  of  the  logarithm  of  the  drag  function1* 

G  =  kBV, 

with  V  expressed  in  meters  per  second ,  and  die 
factor  k  to  take  account  of  differences  in  defini¬ 
tions  and  s'  edards.  Below  600  m/s,  k  =  0.00114; 
but  at  hU  .  r  velocities,  it  is  variable.  G  was 

called  t>  ?.vre  fun-  ion ;  later,  the  syr  Pol  w*s 
changed  G-j  to  denote  that  it  pertains  to  Pro¬ 
jectile  1  ■ 

.n  the  .  oi  3c.-  G,  »  s  argument,  7:’/100, 
v-ruu-  i  0  :c  3  C  (equivalent  to  (817 

m/s  or  ips ) .  ■  .  ,  been  tabulated  for 

V2/100  up  to  25, (XX  ■  ,  .=  m/s  or  5187  fps). 
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The  product  Gv  has  also  been  tabulated  as  a 
function  of  v  in  feet  per  second  up  to  6000  fps.*0 

b.  Projectile  Type  2 

A  sketch  cf  Projectile  T^pe  2  is  shown  in 
Figure  i.-  This  projectile  has  a  6"  boattail  Vi 
caliber  long  and  an  ogivo-conical  head  2.7  cali¬ 
bers  long.  The  experimental  4.7-inch  Shell  IB 
had  this  shape.  However,  the  resistance  firings 
were  actually  conducted  with  two  types  cf  3.3-inch 
shell  :  Type  155  with  a  5  s  boattail  and  T^pe  157 
with  a  7°  boattail.  The  results  were  so  close  that 
they  were  grouped  together  to  represent  one  type. 

These  projectiles  were  fired  from  1922  to  1925 
ct  velocities  horn  655  to  3187  fps.  The  results  cf 
these  firings  were  fitted  by  the  tabulated  drag 
function  G*  (this  is  defined  in  paragraph  3d). 
Tater,  it  was  extrapolated  to  10,000  fps  by  means 
of  theoretical  considerations.  For  the  sake  cf 
greater  smoothness  and  accuracy,  two  revisions 
have  been  made:  G2.i,  tabulated  to  6000  fps,  and 
G2.2,  to  7000  ips.  Ga.;  is  approximately  the  Same 
as  G2,  but  G2.2  is  lower  at  velocities  above  200C 
fps. 

c.  Projectile  Type  8 

A  sketch  cf  Projectile  Type  8  is  shown  in  Fig¬ 
ure  i.  This  projectile  has  a  square  base  and  a 
secant  ogive  2.18  calibers  high.  Its  ogival  radius 
is  10  calibers,  which  is  twice  that  cf  a  tangent 
ogive  cf  tne  sane  height. 

The  drag  function  Gs  was  calculated  from  data 
obtained  by  the  British  Ordnance  Board  from 
tnret  sources: 

(i)  measurements  of  a  model  in  the  National 
Physical  Laboratory  wind  tunnel  at  Mach  numoers 
from  0.22  to  Q.$9. 

'  Resistance  brings  oi  2-pounder  saei-  at 
Mach  numoers  from  0.65  to  2  A 

f  3 >  Resistance  firings  of  3-inch  shell  at  Mach 
numbers  from  0.55  to  4  1. . 

G8  was  first  tabulated  to  5000  ips.  'Rater,  a 
revised  table,  called  Gg.i,  was  extenoed  to  6000 
fps.  The  corresponding  orag  coefficients'  are  ap¬ 
proximately  proDortionai  t~  t-  British  “1940 
Law,”  which  was  Erst  tabulated  to  5000  tps  and 
later  revised  above  3750  fps  and  extended  to  6000 
tps  by  a  power  taw. 

d.  HEAT  Shell,  90-mm,  T108 

A  sketch  of  the  90-mm  High  Explosive  Anti¬ 
tank  Shdl  T108  is  shown  in  Figure  1.  This  pro¬ 
jectile  has  fins  one  caliber  in  span  at  the  end  cf  a 
long  boom  and  a  conical  ogive  about  2  calibers 
high. 

The  drag  cf  this  projectile  was  determined  at 
velocities  from  700  to  2700  fps  by  free  flight  firings 


in  the  large  Spark  Range  of  the  Exterior  Bal¬ 
listic  Laboratory  at  Aberdeen  Proving  Ground, 
The  drag  coefficient  (defined  in  paragraph  3c) 
was  fitted  with  a  series  cf  analytical  expressions, 
which  were  used  in  tabulating  it  as  a  function  of 
Mach  number  from  0.1  to  2.7. 


5.  Trajectories  in  Air 


a.  Equations  oi  Motion 

Under  standard  conditions,  the  equations  cf  mo¬ 
tion  are 


GK  .  ..  '  GH  . 

x  =  ~  — x,  y  =  £-  y  —  g, 


(25) 


where  x  is  the  horizontal  range, 
y — die  altitude, 

H — the  air  density  ratio, 
g — the  gravitational  acceleration  (32.152 
ft/sec2), 

and  dots  denote  derivatives  with  respect  to  time. 
The  initial  conditions  are 


x«  =  0,  y„  —  0,  jc,  =  v0  cos  i>,„  y„  =  v0  sin  i>0, 

where  v0  is  the  initial  velocity, 
v0 — the  angle  of  departure. 

These  equations  cannot  be  solved  analytical'^, 
but  they  can  be  solved  by  a  method  of  numerical 
integration  explained  in  References  21,  22,  23  and 
24. 

h.  Characteristics  of  Trajectories 

(1j  Ballistic  Tobies.  Trajectory  data  have 
been  tabulated  for  Projectile  Types  1  and  2  and 
for  some  bornb  shapes. 

The  “Exterior  Ballistic  Tables  Based  on  Nu¬ 
merical  integration"  (References  7,  25,  and  26; 
pertain  to  Proiectiie  Type  I.  The  meter  is  tne 
unit  of  length.  Tne  trajectories  were  computed 
forward  and  backward  from  the  summi  and  cut 
oft  at  several  altitude.-.  It  C,  is  the  “summitai 
ballistic  coefficient’ ’ct .  trajectoiy  and  yt  the  max¬ 
imum  ordinate  cf  a  partial  trajectory,  the  ballistic 
coefficient  C  of  this  partial  trajectory  satisfies  tne 
relation 


iogioC  =  logicC,  —  0.000,045  y„  (26) 

For  each  y„  C„  and  summital  velocity  v„ 
Volume  I  gives  the  horizontal  range,  feme,  and 
velocity  components  at  the  beginning  and  end 
cf  the  partial  trajectory.  If  the  muzzle  velocitv, 
angle  oi  departure,  arid  ballistic  coefficient  are 
known.  Volume  11  may  be  used  to  find  the 
time  of  flight  to  the  summit,  t»,  the  velocity  at  the 
summit,  the  horizontal  range  to  the  summit,  x*, 
and  the  maximum  ordinate.  Then  C,  can  be  com¬ 
puted  and  trajectory  data  can  be  found  in  Volume 
I  for  the  given  C,  and  v,  at  various  altitudes. 
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However,  terminal  values  can  be  found  more  con¬ 
veniently  from  Volume  III,  which  gives  the  hori¬ 
zontal  range  X.  time  of  flight  T,  angle  of  fall 
id,  and  terminal  velocity  v^,;  the  arguments  are 
the  same  as  in  Volume  II:  the  muzzle  velocity 
varies  from  80  to  920  meters  per  second,  the 
angle  of  departure  from  0"  to  80°  (at  intervals  of 
1°  from  0"  to  .30°,  and  2"  from  30"  to  80°),  the 
ballistic  coefficient  from  1.000  to  17.78  (at  inter¬ 
vals  of  0.05  in  log  C). 

The  "Brief  Exterior  Ballistic  Tables  for  Pro¬ 
jectile  Type  1’’  (Reference  27)  contains  values 
that  were  obtained  by  interpolation  in  Volume 
1 1 1  and  conversion  to  English  units  (yards  and 
feet  per  second) .  The  muzzle  velocity  varies  from 
200  to  3000  feet  per  second  in  100  foot  per  second 
increments.  The  angle  of  departure  varies  from 
15°  to  75"  in  15"  increments.  Obviously,  inter¬ 
polation  for  angle  cf  departure  in  this  table  is 
difficult;  nevertheless,  the  table  is  useful  for  many 


purposes  The  table  of  range  for  an  elevation  of 
45°,  which  is  nearly'  the  maximum  range,  is  repro¬ 
duced  here  (Table  I). 

fhe  "Exterior  Ballistic  Table  for  Projectile 
Type  2"  (Reference  28)  consists  of  three  parts: 
Table  I,  Summital  and  Terminal  Values;  f  able 
IT,  Coordinates  as  Functions  cf  Time  of  Flight; 
and  a  large  graph  of  Maximum  Range  vs.  Ballistic 
Coefficient.  The  tables  have  been  extended  to 
include  data  for  every  15"  from  0"  to  90°  The 
muzzle  velocity  varies  from  1200  to  5200  feet  per 
second,  and  the  ballistic  coefficient  from  0  50  to 
10  00,  The  coordinates  are  given  in  vards,  the 
tinie  in  seconds,  and  the  velocity  in  feet  per  second. 
Table  T  gives  the  time  to  the  summit,  the  hori¬ 
zontal  range  to  the  summit,  the  maximum 
ordinate,  the  time  of  flight,  the  total  range,  the 
angle  cf  fall,  and  the  terminal  velocity.  Table  11 
gives  the  horizontal  range  and  altitude  at  every 
5  seconds  of  time.  Unfortunately,  the  velocity 


TABLE  T 

MAXIMUM  RANGE  FOR  PROJECTILE  TYPE  1 
ELEVATION  45" 

Range  in  yards.  Muzzle  Velocity  (v,)  in  fps 


e, 

Vo 

0.050 

1.122 

0.100 

1.259 

0.150 

1.413 

0.350 

2.239 

0.450 

2.818 

0.550 

3.vv 

0.600 

3.981 

0.660 

4,467 

200 

397 

398 

309 

400 

401 

403 

404 

405 

406 

408 

409 

409 

300 

847 

856 

864 

871 

877 

883 

888 

892 

897 

907 

910 

913 

400 

1419 

1441 

1462 

1481 

1499 

1514 

1528 

1541 

1553 

1.564 

1573 

1582 

1590 

1597 

500 

2077 

2121 

2162 

2201 

2X37 

2270 

2300 

2328 

2354 

2378 

2399 

2418 

2436 

24.52 

000 

2786 

2862 

2934 

3003 

3067 

312ft 

3181 

3232 

3279 

3322 

3361 

3397 

3428 

3461 

700 

3521 

3638 

3749 

3856 

39.58 

4053 

4140 

4222 

4301 

4372 

4438 

4500 

4550 

4606 

800 

4247 

4411 

4724 

4870 

5141 

5267 

5384 

5494 

5595 

5690 

.5777 

5859 

000 

4920 

5136 

5349 

•5556 

5755 

5947 

6133 

6310 

6476 

6634 

6783 

6923 

70.54 

7174 

5489 

5757 

‘H>24 

6286 

6542 

6792 

7037 

7272 

7498 

7714 

7919 

8115 

8299 

8472 

■ffRl 

5925 

6238 

U553 

6866 

7176 

7482 

7784 

8079 

8366 

8644 

8913 

9172 

9420 

9658 

1200 

6251 

6603 

69.58 

7314 

7671 

8027 

8381 

8733 

9077 

9417 

9749 

10073 

10388 

10691 

6521 

7293 

7689 

8087 

8488 

8891 

9292 

9693 

10092 

10489 

10872 

11253 

11626 

6759 

7172 

7594 

8025 

8460 

8905 

9353 

9804 

102.58 

10711 

11165 

11616 

12065 

12508 

1500 

6979 

7419 

7x72 

8336 

8810 

9294 

9787 

10287 

10792 

11301 

11814 

12329 

12845 

13360 

1000 

7185 

7651 

8134 

9141 

9664 

10200 

11302 

11869 

12442 

13021 

1.3605 

14194 

1700 

7381 

7873 

8384 

8912 

9458 

10021 

10599 

11798 

12421 

13056 

13701 

14356 

15019 

1800 

7569 

8088 

8627 

9187 

9767 

10.368 

10989 

11630 

12288 

12967 

13662 

14374 

15099 

15839 

1000 

7752 

8296 

8863 

9454 

10068 

10707 

11371 

12059 

12771 

13505 

14262 

1.5041 

15839 

166.59 

2000 

7931 

8498 

9093 

9715 

10364 

11042 

11749 

12483 

13247 

14037 

14859 

15706 

16580 

17480 

2100 

8106 

8698 

9320 

9972 

10655 

11372 

12121 

12903 

13719 

14.567 

15454 

16372 

17323 

18306 

2200 

8278 

8894 

9.542 

10225 

10942 

11697 

12489 

13319 

14188 

1.5096 

16049 

17038 

18068 

19136 

2:tm 

8447 

9086 

9760 

10473 

11225 

12018 

12853 

13732 

14655 

15624 

16644 

17706 

18820 

19975 

24'  0 

8612 

9274 

9975 

10718 

11.505 

12336 

1321s 

14143 

15121 

16152 

17239 

18375 

19574 

20821 

2500 

8774 

9459 

10188 

10961 

11732 

12652 

13575 

14552 

15585 

16679 

17834 

19048 

20)130 

21675 

2600 

i  8935 

9642 

10397 

11201 

12056 

12965 

13931 

14958 

16048 

17205 

18430 

19725 

21095 

22539 

2700 

KSZI 

9823 

10603 

11438 

13237 

13275 

14285 

15362 

19027 

20405 

21866 

23412 

2800 

9247 

10001 

10807 

11672 

12595 

13582 

14637 

15764 

16970 

18257 

19625 

21091 

22645 

24294 

2000 

9399 

10177 

11010 

11903 

12861 

13886 

14987 

16165 

17429 

18782 

20225 

21778 

23428 

25184 

3000 

9549 

10351 

12132 

13123 

14187 

15334 

16563 

r  .385 

19303 

20823 

22464 

24213 
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TABLE  I  (Continued) 

MAXIMUM  RANGE  FOR  PROJECTILE  TYPE  1 
ELEVATION  45° 


Range  in  yards.  Muzzle  Velocity'  (v0)  in  fps 


Log  Ci.. 

c,.. 

Vo 

0.650 

4.407 

0.700 

6.012 

0.750 

5.623 

0.800 

6.310 

0.850 

7.079 

0.900 

7.943 

0.950 

8.913 

1.000 

10.00 

I. 050 

II. 22 

1.100 

12.59 

m 

1.200 

15.85 

1.250 

17.78 

200 

409 

409 

409 

410 

410 

410 

411 

411 

412 

412 

412 

413 

413 

413 

300 

910 

913 

915 

917 

918 

919 

921 

923 

924 

925 

925 

926 

927 

928 

400 

1590 

1597 

1604 

1611 

1&4 

1619 

1624 

1627 

ELLjI 

1634 

1637 

1639 

1641 

1643 

500 

2436 

2452 

2467 

2480 

2491 

2501 

9511 

2519 

2528 

2533 

2540 

2545 

2550 

2555 

600 

3428 

3461 

3489 

3513 

3535 

3555 

3575 

3591 

3605 

3619 

3631 

3649 

3659 

700 

4556 

4606 

4651 

4695 

4733 

4766 

4800 

4829 

4855 

4879 

4899 

4918 

4934 

4951 

Bltil 

5777 

5859 

5934 

6001 

6064 

6120 

6172 

6219 

6260 

6298 

6333 

6364 

6394 

/TilO 

i 

900 

7054 

7174 

7284 

7389 

7484 

7573 

7652 

7725 

7791 

7852 

7907 

7958 

8003 

044 

1000 

8299 

8472 

8634 

8786 

8928 

9060 

9180 

9292 

9394 

9489 

9575 

WT87 

1100 

9420 

9658 

9882 

10094 

10295 

10484 

10660 

10825 

10979 

11120 

11252 

11373 

11485 

11588 

1200 

10388 

10691 

10984 

11286 

11538 

11794 

12038 

12288 

12487 

12691 

12882 

13062 

13228 

13382 

1300 

11253 

11628 

11990 

12344 

12688 

13020 

13339 

13644 

13936 

14213 

14476 

14723 

14958 

15175 

1400 

12066 

12508 

12945 

13375 

13796 

14207 

14604 

14991 

15364 

15722 

16067 

16394 

16097 

1500 

12845 

13360 

13872 

14381 

14881 

15374 

15857 

16331 

16790 

17235 

17867 

18474 

18852 

1600 

|R?S!S 

14194 

14783 

15370 

15956 

16536 

17109 

17674 

18227 

18766 

19291 

19798 

ra 

1700 

14356 

15019 

15686 

16357 

17030 

17702 

18371 

19031 

19682 

20324 

20949 

21559 

1800 

15839 

16588 

17344 

18108 

18875 

19643 

20406 

21102 

21913 

1900 

15839 

16659 

17491 

18338 

19198 

20063 

20936 

21809 

22678 

23543 

2000 

16580 

18399 

19340 

20299 

21270 

22253 

23243 

24222 

25222 

2100 

17323 

18306 

19315 

20353 

21415 

22497 

23598 

24710 

25829 

m5  7 

2200 

18068 

19136 

21380 

22551 

23748 

24972 

21S217 

2300 

18820 

19975 

21177 

22421 

23707 

25027 

26385 

27771 

2400 

19574 

20821 

22124 

23478 

24883 

26334 

27832 

29372 

2500 

£ 

21675 

24553 

26083 

27673 

29321 

31028 

2600 

21095 

22539 

24056 

29048 

30861 

32748 

2700 

21866 

23412 

25045 

30458 

32446 

|  2800 

22645 

24294 

28040 

27890 

29849 

31910 

34091 

j  2900 

23428 

25184 

27049 

29037 

31183 

33410 

35797 

|  3000 

24213 

26083 

28080 

34944 

37558 

_ s 

TABLE  II 

MAXIMUM  RANGE  OF  PROJECTILETYPE  2 
Range  in  Units  of  10  Yards 


.  Log  C» 

Wlrwntv  I - 1  .  i - : - 


Mu* He  velocity 
fp. 

-0.3 

-0.2 

-0.1 

0.0 

0.2 

0.4 

o.e 

0.8 

1.0  j 

1200 . 

709 

861 

1004 

1128 

1227 

1308  i 

1600 . 

851 

1088 

1335 

1584 

€821 

2000 . 

1001 

1322 

1990 

2091 

2499 

2400 . 

1143 

1561 

2071 

2659 

3295 

3983 

2800 . 

759 

905 

1078 

1284 

1806 

2479 

4322 

5359 

3200 . 

819 

983 

1182 

1426 

2062 

2926 

4074 

5589 

7006 

3600 . 

877 

1061 

1287 

1571 

2335 

3438 

5061 

7201 

9034 

4000 . 

934 

1139 

1394 

1723 

2629 

4066 

6329 

9106 

11356 

4400 . 

992 

1218 

1504 

1874 

2950 

4834 

7917 

11335 

14019 

4800 . 

1617 

3342 

5845 

9836 

13888 

17015 

5200 . 

1735 

2220 

3795 

7196 

12091 

16752 

20333 

components  are  not  tabulated.  The  graph  has  a 
curve  of  maximum  range  vs.  ballistic  coefficient 
for  every  400  feet  per  second  muzzle  velocity;  it 
also  lists  the  angle  of  departure  for  maximum 
range.  A  table  the  maximum  nuige  is  included 


in  this  pamphlet  (Table  11). 

Some  trajectories  have  been  computed  for  Pro¬ 
jectile  T^pe  8  at  an  angle  <£  departure  of  45°. 
The  range,  which  is  nearly  the  maximum  range, 
is  tabulated  here  (Table  111). 
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TABLE  III 

MAXIMUM  RANGE  FOR  PROJECTILE  TYPE  8 
ELEVATION  45” 

Range  in  Units  of  10  Yards 


Mu*ile  Velocity  »p» 

-0.2 

0.0 

0.2 

0.4 

0.6 

0.8 

1.0 

1200  . 

521 

665 

815 

961 

1091 

1201 

1288 

2000 . 

680 

923 

1223 

1581 

1934 

2403 

2826 

2800. 

822 

1168 

1650 

2287 

3124 

4231 

5251 

3600  .  . 

964 

1432 

2145 

3162 

4865 

6936 

8741 

4400. 

1108 

1720 

2726 

4470 

7492 

10884 

13661 

5200  . 

1258 

2040 

3424 

6604 

11693 

16582 

20319 

When  bombs  were  dropped  at  low  speeds,  tra¬ 
jectory  data  based  on  the  Givre  drag  function  Gi 
were  accurate  enough  for  calibrating  the  sights. 
In  fact,  the  descending  branches  of  the  trajectories 
in  Volume  I  of  the  “Exterior  Ballistic  Tables” 
could  be  used  for  this  purpose.  Actually,  some 
“Bomb  Ballistic  Reduction  Tables”  were  pub¬ 
lished,  which  were  more  convenient  for  this  pur¬ 
pose.  However,  at  transonic  speeds,  the  drag 
functions  of  modern  bombs  differ  considerably 
from  Gt.  Drag  coefficients  have  been  determined 
for  several  shapes,  but  three  cf  them  have  been 
selected  as  bases  of  “Bomb  Ballistic  Tables”  (Ref¬ 
erence  10):  drag  coefficient  A  was  determined 
from  wind  tunnel  tests  cf  a  shape  designed  by  the 
National  Advisory  Committee  for  Aeronautics,  and 
is  probably  the  lowest  drag  coefficient  that  any 
bomb  will  have  ;  drag  coefficient  B  was  determined 
from  range  bombing  of  the  3000-pound  Demoli¬ 
tion  Bomb  Ml  19,  and  is  suitable  for  other  mem¬ 
bers  of  the  series  including  the  750-pound  Bomb 
Ml  18  and  the  10,000-pound  Bomb  M121;  drag 
coefficient  C  was  determined  from  resistance  fir¬ 
ings  of  the  75 -mm  Slug  Mark  I,  and  is  probably 
the  highest  coefficient  that  any  bomb  will  have. 

The  arguments  of  the  “Bomb  Ballistic  Tables” 
are  release  altitude  ( 1000  to  80,000  feet  above  sea 
level),  true  air  speed  cf  release  in  level  flight  (250 
to  2000  knots),  and  performance  parameter  {0 
to  1.5).  The  performance  parameter  P  is  inverse¬ 
ly  proportional  to  the  ballistic  coefficient;  it  is  de¬ 
fined  by  the  formula 

P=100d7m,  (27) 

where  d  is  the  maximum  diameter  df  the  barb  in 
feet. 


m — the  mass  of  the  bomb  in  pounds. 

The  tables  give  the  horizontal  range  from  release 
to  impact  (feet),  the  time  from  release  to  impact 
(seconds),  the  striking  velocity  (ft/sec),  and  the 
striking  angle  (degrees).  The  drag  coefficients 
are  tabulated  as  functions  cf  the  Mach  number  in 
the  appendix;  they  are  graphed  in  the  intro¬ 
duction. 

(2)  Graphs.  Figures  3,  4  and  5  are  graphs  of 
the  total  time  of  flight  for  an  elevation  of  45”.  The 
time  cf  flight  is  useful  in  estimating  wind  effects. 

For  antiaircraft  firing  and  for  firing  to  the 
ground  at  different  altitudes,  points  along  the  tra¬ 
jectory  are  needed.  Figures  6  and  7  show  the 
trajectories  computed  with  G2  and  GH  at  angles  of 
departure  of  45”  and  60°  (there  is  a  different 
graph  for  each  muzzle  velocity,  each  elevation, 
and  each  Projectile  Type).  Trajectories  based 
on  Gi  are  not  included  because  no  antiaircraft  pro- 
jectiies  have  a  blunt  head.  The  T108  HEAT 
Shell  is  usually  used  only  in  direct  fire  at  short 
ranges.  Trajectory  data  for  other  elevations  can 
be  estimated  from  those  at  45°  or  60°  by  assuming 
that  the  slant  range  is  the  same  for  a  given  time  cf 
flight. 

Unlike  the  vacuum  trajectories,  in  air,  the 
horizontal  range  to  the  summit  is  more  than  hall 
the  total  range,  and  the  time  to  the  summit  is  legs 
than  half  the  total  time;  however,  the  maximum 
ordinate  is  approximately  (g/8)T2.  In  air,  the 
terminal  velocity  is  less  than  the  initial  velocity 
and  the  angle  cf  fall  is  greater  than  the  angle  cf 
departure.  At  a  given  elevation,  the  total  ryige 
increases  less  rapidly  than  the  square  of  th* 
muzzle  velocity;  at  a  given  muzzle  velocity,  the 
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range  is  a  maximum  at  ail  angle  somewhat  differ¬ 
ent  from  45°  :  lower  for  small  ballistic  coefficients, 
higher  for  large  ones.  However,  the  range  at  45° 
is  nearly  as  much  as  the  maximum. 

6.  Siacci  Tables 
a.  Spin-stabilized  Projectiles 

(1)  Tabular  Functions.  In  1880,  the  Italian 
Col.  Siacci  introduced  the  Space,  Time,  Altitude, 
and  Inclination  Functions,  (Reference  29) ,  which 
simplify  the  calculation  of  trajectoiy  data  under 
certain  restrictions.  In  particular,  their  use  will 
be  explained  for  three  conditions:  ground  fire  at 
low  superelevations,  antiaircraft  fire  (where  the 
superelevation  is  small),  and  aircraft  fire.  These 
and  some  auxiliary  functions  based  on  Gi,  G2.S 
(the  second  revision  of  G2)  and  G8.i  (the  first 
revision  of  G8)  have  been  tabulated  for  the  Air 
Force  (References  30,  31  and  32).  Since  these 
Siacci  tables  were  arranged  primarily  for  use  in 
computing  aircraft  firing  tables,  they  were  called 
Aircraft  Tables;  however,  they  can  be  used  for 
other  pmposes,  as  explained  below,  in  lieu  of  other 
fonns.  In  these  Aircraft  Tables  the  argument  is 
the  component  u  along  the  line  of  departure  of  the 
velocity  relative  to  the  air,  which  is  assumed  mo¬ 
tionless.  The  following  functions  and  some  of 
their  derivatives  with  respect  to  u  are  tabulated: 


Space 

p-fU 

du 

(28) 

Ju 

G(u)  ’ 

Time 

r+ 

II 

c 

du 

(29) 

uG(u)  ’ 

Inclination  * 

ru 

*=J 

•Ml 

g  du 

(30) 

u*G(u) 

Altitude"  q=J'  "cfyuj — ’  (31) 

Velocity  increment  Alt  =  G(u),  (32) 

(Kdj  =  16.4 per  rad8), 

Damping  coefficient  c"  =  G(u)/2u,  (33) 

The  upper  limit  of  integration,  U,  is  6000  fps  for 
Projectile  Type  1  and  7000  fps  for  Types  2  and 
8.  Kds  is  the  yaw-drag  coefficient,  defined  by 
the  formula 

KD  =  KD#(l-fKD  82),  (34) 

where  Kd*  is  the  drag  coefficient  for  0  yaw, 

4 — is  the  yaw  in  radians. 

(2)  Ground  Fire  at  Low  Elevations 

(a)  Formulas.  Although  the  Aircraft  Tables 
are  not  accurate  in  computing  complete  trajectories 
at  high  elevations,  they  can  be  used  to  obtain  tra¬ 
jectory  data  at  superelevations  below  15"  (for  the 
derivation  and  application  of  the  following  for¬ 
mulas,  see  Reference  33).  Here,  we  shall  take 

jj/po — as  the  ratio  cf  air  density  to  normal,  as¬ 
sumed  constant  for  each  problem, 

a/a* — the  ration  of  the  velocity  of  sound  to 
normal,  assumed  constant  for  each 
problem. 

F  F  is  the  temperature  in  degrees  Fahrenheit, 

a/a*  =  V(459-f-F)/518.  (35) 

Then  the  horizontal  range  is 


^  =  (Cp0/p)  cos  $0  [p(ua*/a)  —  p(u0a0/a)],  (36) 

the  time  of  flight  is 

tf=  (CpoSo/pa)  [t(ua*/a)  —  t(u*a*/a)],  (37) 

the  altitude  is 

y  =  xtan30  — (Cp^/pa)*  [q(ua*/a)  —  q(u*a*/a)]  -f  x(Cp*a*/pa)  sec$0  fq1(uia<i/a)],  (38) 
and  the  slope  of  the  trajectory  is 

tan  $  =  tan30— (Cp*a*VPa*)  ^  $0  [qt  (ua*/a)  —  qi  (u*a*/a)].  (39) 


The  angle  of  departure  for  y  —  0,  may  be  found  by  formula 


sin  i3. 


Cp*a*  r q(ua*/a)  —  q(u*a*/a) 
?a  Lp(ua*/a)  —  p(u*a*/a) 


qi(u*a*/a)J. 


If  p/p0  —  1  and  a/a*  =  I,  let  p  =p(u),  p0  =p  (u*),  etc. 


(40) 


•In  most  places,  the  altitude  and  inclination  functions  are  defined  as  twice  these  integrals. 
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Then  these  formulas  can  be  expressed  more 


simply  : 

x  =  Ccos9,  (p  —  p«),  (36a) 

tf=C(t—  to),  (37a) 

y  =  x  tan  $0  —  C2(q  — q0) 

CsecS0xq,0)  (38a) 

tan  S  —  tan  9,  —  C  sec  9,  (qi  — qio).  (39a) 
and,  for  y  ~  0, 

sin  ^  C  —  —  qi0  1  (40a) 

Lp  Po  J 


(h)  Example 

Data:  Tank  Gun, 
M339, 

Caliber 
Mass 

Muzzle  Velocity 
Form  Factor 
Air  Density  Ratio 
Sound  Velocity  Ratio 


M1A2 ;  Shot,  AP. 

3.000  in. 

14.58  lb. 

3200  fps 
1.164  on  Ga.i 
1.000 
1.000 


76-mm, 


Problem:  To  find  the  range,  time  of  flight,  angle 
of  departure,  and  angle  of  fall  where 
the  velocity  is  2000  fps. 


Solution  :  From  the  Aircraft  Table  Based  on  Gg.i : 


u 

P 

t 

q 

<h 

2000 

23079.3 

6.5132 

489.0 

.06523 

3200 

15425.1 

3.4997 

161.6 

.02637 

Diff. 

7654.2 

3.0135 

327.4 

.03886 

CSJ  =  14.58/1.164 

(3.000) 2  = 

1.392  (Eq.  22) 

sin  9,  = 

1.392  (327.4/7654,2  - 

-.02637) 

— 

.02283  (Eq.  40a) 

9.  =  23.26  mils,  cos  9,  =  .99974,  sec,  i$0  = 
1.00027 


tan  ■$„  =  .02284 

1.392(1.00027)  (.03886)  =  .0541  1 

tan  9  =  — .03127  (Eq.  39a) 

9  =  —31.84  mils 

x  —  1  392 (.99974)  (7654  2)  = 

10652  ft  (Eq.  36a) 

t,  =  1.392  (3.0135)  =  4.19  sec  (Eq.  37a) 


Note:  Since  u  is  the  component  cf  the  velocity 
along  the  line  cf  fire,  the  tangential  velocity  is 
v  =  u  cos  S0/cos  9. 

In  this  example,  v  =  2000  (.99974 )/.99951  = 

2000.5  fps. 

(3)  Antiaircraft  Fire 

(a)  Formulas.  These  tables  can  also  be 
used  for  the  flat  part  of  antiaircraft  trajectories, 
where  the  superelevation  is  less  than  15°  (for  the 
derivation  and  application  cf  the  following  formu¬ 
las,  see  Reference  33).  Let  p/p0  and  possibly  C 
be  known  variable  functions  of  x,  but  assume  a 
constant  value  cf  a/a*  (either  1  or  an  average). 
Since  p/p0  is  a  function  cf  y,  denoted  by  H,  its 
relation  to  x  can  be  found  by  the  approximate  sub¬ 
stitution 

y  «=«  x  tan -Sc,  (41) 

Let  m  =  tan  9,  m,  =  tan  9,. 

Then  the  following  trajectory  data  can  be  found 
by  Simpson's  rule  or  some  other  method  of  nu¬ 
merical  integration  : 


x 

p  (ua<,/a)  =p(u0a„/a)  +  sec^J'^  dx,  (42) 

o 


X 

rdx 

—  .  (43) 

o 


m 


f  gsec 

=  m°~  J  — 


sec2  9, 


dx, 


(44) 


(45) 


Check:  xtan  •So  —  10652  (.02284)  =  243.29 

—  <3.1  (q  —  q0)  =  (1-392)*  (327.4)  =  —  634.40 

Cgl  sec  ■&„  X  qio  =  1.392  (1.00027)  ( 10652)  (.02637)  =  39111 

y=  0.00  (Eq.38a) 
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(b)  Example 


Ditto:  Automatic  (IVllGun.  40-mm,  M2A1;  Shot,  AP,  MSI 

Muxile  Velocity 

2870  fpe 

Balliatic  Coefficient 

Q1574  on  G«.t 

Angle  cf  Departure 

:  1200  mil* 

Air  Denaity  Ratio 

:  H-e^r 

Sound  Velocity  fetio 

1,000 

Problem:  To  find  the  time  of  flight  rad  altitude  for  horizontal  ranges  from  0  to  2400  ft. 

Soluticn:  (See  remark*  below  aid  uae  the  Aircraft  Table  Bleed  on  G,.,) 

H 

J~Hdx 

p(u) 

U 

‘f 

/2813.iy 

r/ 2613.1  \*Hx 

m 

H 

m 

u 

V  u  ) 

n  u  ) 

tt 

fpt 

aec 

H 

H 

3 

4 

6 

6 

7 

8 

9 

10 

ii 

12 

0 

300 

600 

.00000 

.00903 

.01987 

1.0000 

.9774 

.9653 

0.00 

17439.6 

2870.0 

0.9105 

0.000 

0.8290 

0 

2.4142 

0 

588.49 

20109.5 

2446.4 

1.0681 

0.594 

1.1408 

591 

2.3952 

1443 

900 

1200 

.02980 

.03973 

.9337 

.9126 

1146.76 

22660.1 

2060.6 

1.2681 

1.295 

1.6081 

i.4i6 

2.3687 

2872 

1500 

1800 

.04966 

.05960 

.8919 

.8718 

1681.96 

25096.6 

1716.3 

1.5225 

2.132 

2.3180 

2593 

2.3306 

4282 

2100 

2400 

.06953 

.07946 

.8521 

.8328 

2193.26 

27424.2 

1413.7 

1.8484 

3.143 

3.4166 

4314 

2.2755 

5664 

Remarks  : 

Col.  1.  The  interval  of  x  should  be  chosen  small 
enough  so  that  the  integrations  may  be  performed 
accurately  by  Simpson’s  rule:  not  more  than  600 
feet. 

Col.  2.  M  =  log,,  e=  0.43429.  h  =  0.000,031,- 
58  ft-1.  The  estimate  of  y  =  x  tan  $<>•  A 
slightly  smaller  factor  would  be  better  at  the 
longer  ranges. 

Col.  3.  For  convenience,  the  formula  H  =  anti- 
cologio  Mhy  was  used.  The  density  function 
H(y)  has  been  tabulated  (Reference  11)  ;  it  gives 
the  same  values.  The  ICAO  “Standard  Atmos¬ 
phere — Tables  and  Data"1  gives  slightly  different 
values. 

Col.  4.  Simpson’s  rule  was  used,  but  the  trap¬ 
ezoid  rule  gives  2193.30  for  the  last  value. 

Col.  5.  The  first  value  was  found  in  the  Aircraft 
Table  for  Uo  =  2870  fps.  The  other  values  were 
computed  by  Eq.  (42)  with 

sec  Vc8.i  =  2.6131/0.574  =  4.5524. 

Col.  6.  The  values  of  u  after  the  first  were  found 
in  the  Aircraft  Table  by  inverse  interpolation, 
using  the  first  two  terms  cb  Taylor’s  series. 

Col.  7.  sec  =2.6131.  This  is  multiplied  by 
1000  to  avoid  writing  zeros. 


Col.  8.  Eq.  (43) .  Simpson’s  rule  gives  3.131 
for  the  last  value,  but  the  trapezoid  rule  was  used 
in  order  to  obtain  intermediate  values. 

Col.  9.  The  square  of  Col.  7. 

Col.  10.  Note  the  factor  10®  in  these  values. 
Simpson’s  rule  gives  4259  for  the  last  value,  but 
the  trapezoid  rule  was  used  so  as  to  obtain  more 
values. 

Col.  11.  Computed  by  Eq.  (44)  with 

m,  =  tan  ■So,  g  =  32.152ft/sec*. 

Col.  12.  Computed  by  Eq.  (45).  Simpson’s  rule 
gives  5666  ft  for  the  final  value. 

(4)  Aircraft  Fire 

(3)  Formulas.  These  tables  are  especially 
useful  for  aircraft  fire  (for  ths  derivation  and  ap¬ 
plication  of  the  following  formuas,  see  References 
34  and  35) .  Before  using  them,  the  initial  con¬ 
ditions  must  be  determined.  I^et 
A-tee  the  azimuth,  measured  from  the  direc¬ 
tion  of  motion  of  the  airplane, 

2 — the  zenith  angle,  measured  from  the  upward 
vertical, 

v„ — the  muzzle  velocity  (relative  to  the  gun), 
w — die  true  air  speed, 

Uo — the  initial  velocity  relative  to  the  air, 

$0 — the  initial  yaw  relative  to  the  motion  cf  the 
air. 
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The  initial  velocity  and  yaw  may  be  found  by  the 


formulas 

u,*  =  v.*  -f-  w*  -4*  2wv0  sin  Z  cos  A,  (46) 

8,*  =  (1  —  sin*  Zcos*  A)  w*/uo*  (rad)1.  (47) 
In  forward  tiring,  A  =  0;  Z  =  1600  mils, 
u,=v0-t‘w,  (46a) 

5.  =  0.  (47a) 

In  rearward  filing,  A  =  3200  mils,  Z  =  1600 
mils, 

u,  =  v„  —  w,  (46b) 

8o=0.  (^7b) 

The  air  density  raho  H  and  the  sound  velocity 


ratio  a/ao  may  be  taken  as  constants.  Approxi¬ 
mately,  if  Y  is  the  altitude  of  the  airplane  above 
sea  level, 

H  -  h  =  0.000,031,58  ft-*,  (48) 

a/a.  =  e-«,  k  =  0.000,003,01  ft-*.  (49) 

If  s,  is  the  standard  stability'  factor,  which  cor¬ 
responds  to  the  muzzle  velocity,  muzzle  spin,  and 
standard  surface  air  density,  the  initial  stability 
factor  is 

s»  =  s,v,*/Uo*  H.  (SO) 


If  h  is  the  yawing  moment  damping  factor, 
x— the  cross  wind  force  damping  factor, 

c’=  (h  +x)/2u,  (51) 

the  loss  in  velocity  due  to  yaw  is 
a..  _  *«*a  ( s«  — 1/2)  A'  u  (uo/a) 

AU  “  "IT  (s.— l)Cc'  +  c"  (uoa./a)  (52) 
Then  the  effective  initial  velocity  is 


u,  =  u,  —  Au.  (S3) 

Now  we  can  find  the  distance  measured  along 
the  line  db  departure  to  a  pint  directly  above  the.* 
projectile, 

P  =  (C/H)  (p(uao/a)  —  p  (iMo/a)],  (54) 

the  drop  cf  the  projectile, 

Q—  (Cao/Ha )*  [q(uao/a)  —  q(u«a./a)] 

—  (Cao/Ha)  Pqj  (ua/s),  (55) 

and  the  time  of  flight  is, 
tf  =  (Cao/Ha)  {t(uao/a)  — t  (uA/a)].  (56) 
The  horizontal  and  vertical  components,  if  de¬ 
sired,  may  be  found  by  the  formulas 

x  =  Pcos6„,  (57) 

y  =  P  sin  —  Q,  (58) 

with  the  angle  of  departure  determined  by  the 
formula 


t  4  «  _ _ fos'Z _ 

'  (sin  Z  sin  A)="‘  (sinZ  cos  A -f  w/vc)* 

or 

cot*  6.=  (tan  Z  sin  A)*  -|-  (tan  Z  cos  A  -j-  w/v0cosZ)*. 


(59) 

(59a) 


Illustrative  example  is  given  on  the  following  page. 
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(b)  Example 
Data: 

Azimuth 
Zenith  Angle 

Muzzle  Velocity 
Air  Speed 

Altitude 

Standard  Stability  ESdrr 
Damping  Factor*  h  +  x 

Ballistic  Coefficient 


1200  mils 
800  mils 
3000  fps 
1000  fps 
30,000  ft 
3.66 

10.6  sec-1 
at  3000  fps 
0.235  on  Gs.i 


Problem:  To  find  the  time  cf  flight  and  drop  for 
a  Siacci  range  (P)  at  1200  ft 

Solution:  (Use  the  Aircraft  Table  Basal  on 
cos  A  ss  0,38268 


sin  Z  =  0.7JZL1 
sin  Z  cos  A  =  0.27060 


2wv„  =  2  x  1000  x  3000  =  6,000,000 
v,1  =  (3000)* 


w*  =  (1000) 3 

2wv0  sinZ  cos  A 

14* 

Uo 


=  9,000,000 
=  1,000,000 
=  1,623,600 


=  11,623,600  (Eq.46) 
=  3409.3  fps 
80*  =  0.72940  x  1,000,000/11,623,600  = 
0.06275  (Eq.  47) 

3s,  =0.2505  rad  (2552mils),  (1  rad  = 
lQ18.6mils). 

hY  =0.000,031,58  X  30,000  =  0.9474 

H  =  0.3877  (Eq.  48) 

kY  =  0.000,003,01  x  30,000  =  0.0903 


a/a.  =  0.SO7(Eq.  49) 
s.  =  3.66  x  9,000,000/1 1,623,600  x  0.3877  = 
7.3L  (Eq.  50) 

c'  =  10.6/2  x  3000  =  0.001.77ft-1 
u^/a  =  3409.3/0.9137  =  3731.3  fps 

A'u(us>ao/a)  =  1.4003  sec”1 

(from  Aircraft  Table) 


c"(i4a«/a)  =  2288  x  10"*  ft-1 
(from  Aircraft  Table) 


u.  =  340921  -  206.5  =  3202.8  fps  (Eq.  53) 
Cm/H  =  0235/0.3877  -  0.6061 
CriBa/H*  =  0.6061/0.9137  =  0.6633 
(Cm4/H»)‘  =  (0.6633)*  =  0.4400 
PH/Cn  =  1200/0.6061  =  1979.9  ft 
uA/a  =  3202.8/0.9137  =  3505.3  fps 
Find  p(v4S4/a)  in  Table,  then  p(ua*/a)  by  Eq. 
54. 


ua./a 

P 

t 

q 

9* 

3174.3 

15580.5 

3.548 

165.7 

.02686 

3505.3 

13600.6 

2.9$5 

HS.4 

.02114 

1979.9 

0.598 

47.3 

.00572 

t,  =  0.6633  x  0.933  =  0.3S3sec  (Eq.  56) 

0.4400  x  47.3  =  20.8 

0.6633  x  1200  x  0.02114  =  16.8 

Q  =  4.0  ft  (Eq.55) 

(c)  Interpolation,  The  derivatives  listed  in 
the  Aircraft  Tables  should  be  used  for  inter- 
pdLatiiig  backward  or  forward  fern  the  nearest 
tabular  value  of  the  argument.  Er  instance,  in 
the  last  example,  it  is  desired  to  find  p  for  u  — 
3505.3  and  then  to  find  u  for  p  =  15580.5.  Errst, 
the  table  gives 


at 


u=  3510, 

dp  _ 

du 


=  -  5.£3 


p  =  13572.81, 


Hence,  Au  =  —  47,  Ap  =  -f  27.84, 

and,  at  u  =  3505.3,  p  =  13600.65. 

Then  we  have  p  =  15580.5 
and  the  table  gives 


p  =  15606.2,  ^-=  _  6.044  et  u  =  3170. 

du 

Hence,  Ap  =  —  25.7  Au  =  43, 


and 


u  =  3L74.3. 


=  0.06275  nqn7 _ 6.81x1.4003 _ 

1  XU  6.31  x  0.235  x  0.00177  +  0.000,022,88 

=  005733  0.002,625a+ 0.000.023  =  2055  ^  <^S2> 
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b.  Fin-stabilized  Projectiles 

(1)  Tubular  Functions.  Reference  6  contains 
a  table  of  modified  Siaeci  functions  based  on  the 
drag  coefficient  for  the  90-mm  High  Explosive 
Antitank  Shell  T108,  which  has  rigid  fins.  These 
tables  can  be  used  for  other  similar  projectiles, 
although  they  are  not  suitable  for  all  fin-stabilized 
projectiles.  The  following  functions  and  their 
derivatives  are  tabulated  as  functions  of  Mach 
number  M  : 


2.7 


Space 


Time 


S-  f— d-M. 

J  MKd(] 


M 

2.7 


(M)  ’ 


T 


=  f 


6M 


M 

2.7 


M2Kd(M)’ 


Inclination  I  = 


g_r _ dM 

a J  M*Kd 


M 

2.7 


Altitude  A 


-  CMM) 

-  J  MKd 
M 


(M)’ 

I(M)dM 
(M)  • 


(60) 


(61) 


(62) 


(63) 


The  factor  in  the  inclination  function  is 

g/s Lo  =  32.154/1120.27  =  0.0287  sec"1. 


y  =  x  tan  9, 


_  aoPo2  sec2  40 


a.2p2f2 


[A(M)  —  A(M0)J 


tan  4  =  tan  40  -  (I(M)  ~I(M0)]. 

The  angle  of  departure  for  y  =  0  may  be  found  by 


sin  240 


_  a.,Po  [~*  A(M)  —  A(M0) 
a2PTL  S(M)  -  S(M0) 


(2)  Application 

(a)  Formulas.  The  following  formulas  are 
applicable  for  computing  trajectory  data  at  super¬ 
elevations  below  15°. 

Let 

m  be  the  mass  of  the  projectile  in  pounds, 
i — the  form  factor  of  the  projectile, 
d — the  maximum  diameter  of  the  projectile  in 
feet, 

v — the  velocity  of  the  projectile  in  feet  per 
second, 

v0 — the  muzzle  velocity  in  feet  per  second, 
a — the  velocity  of  sound  in  air  in  feet  per  second, 
j> — the  density  cf  air  in  pounds  per  cubic  foot, 
9 — the  angle  of  inclination  cf  the  trajectory, 

•$<> — the  angle  of  departure. 

Also  let 

a,=  1120.27  fps, 
po  =0.07513  lb /ft8, 

M0  =  vo/a,  (64) 

M  =  v  cos  4/a  cos  (65) 

y  =  j>0id2/m.  (66) 

Then  the  horizontal  range  in  feet  is 

x  =(po/pt)[S(M)  — S(M.)],  (67) 

the  time  cf  flight  in  seconds  is 

tf=  (p.sec  V«PT)  JT(M)  -  T(M.)],  (68) 
the  altitude  in  feet  is 


+  x 


aoPo  sec2  40 
a2pr 


I(Me). 


(69) 


(70) 

the  formula 


(71) 


If  p  =  po  and  a  =  ao,  let  S  =  S(M),  S0  = 

S(M0),  etc.  Then  these  formulas  can  be  ex¬ 
pressed  more  simply: 


x  =  (1/y)(S  —  S„),  (67a) 

tf  =  (sec.40/aoY)(T  -T.),  (68a) 


y  =  x  tan  40 


sec* 40  ,  A  A  ,  , 

- (A  —  A<>)  +x 


see2  4„ 
aoY  ° 
(69a) 


tan  4  =tan  40  -  (I  —  I,),  (70a) 

and  for  y  =  0, 


<7I*> 
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(b)  Example 


Data:  Gun,  90-mm;  Shell,  HEAT ,  T108E20 


Caliber 

Mass 

Muzzle  Vdocity 
Form  Factor 

Air  Density 
Sound  Velocity 


:  029525  ft 

:  142  lb 
:  2400  fps 
:  1.00  on  Gnoc 

:  0.07513  lb/ft4 
:  112027  fps 


Problem :  To  find  the  velocity,  time  cf  flight,  angle 
of  departure,  and  angle  cf  fill  at  a  range 

cf  3000  feet 


Solution :  From  the  Table  cf  the  Siacd  Functions 
based  on  Shell,  90-mm,  HEAT,  T108, 

MS  T  A  I 

1.73650  3.091%  1.427502  .02537  .01923 

2.14234*  170836  0.709983  .00670  .00851 

EiE  08360**  0.717519  .01867  .01072 


4  =  2400/112027  =  2.14234  (Eq.  64) 

y  =  0.07513(  1.00)  (029525)*/142  = 
0.000,4612  ft-1  (Eq.66) 

Ut  =  112027(0.000,4612)  =  0.51667  sec'1 

a.-r1  =  0.51667(0.000,4612)  = 

0.00023829  ft-1  sec-1 

**S  -S.  =rx  =0.000,4612(3000)  = 

1.38360  (Eq.  67a) 

sin  241,  s  ftf  &67  -  0.00851  J 

=  0.01928  (Eq.  71a) 

241 5_=  19.64 mils,  $0  =  9.82  mils,  cos  4),  = 
0.99995, 

sec-S.  =  1.00005,  sec1  a,  =  1.00010 
tan  •S,  =  0.00964 

-  (0.01072)  =  -0,02075 

tan  #  =  -0.01111  (Eq,  70a) 

$  =  —  1 1.32  mils,  cos  4)  =  0.9994 

v=  112027  (173650)  (0.99995)/(0.95994) 

=  1945.37 fps  (4.65) 

1  00005 

U~  -^^  (0717519)  =  1.389  sec  (Eq.68a) 


Check:  x  tan «.  =  3000  (.00964)  =28.92 ft 
sec*  *.(A  —  A.)  1.00010(0.01867) 

a.y*  ~  0.00023829 

=  -78.36  ft 

x  sec* T  3000(1.00010)  , 

0.51667 - (0-°°851) 

=  +49.42  ft 

y  =  28.92  +  49.42  -  78.36  =  —0.02  ft 
(Eq.  69b). 

7.  Diffexential  Effects 
Some  of  the  effects  on  range  and  deflection 
caused  by  variations  from  standard  conditions 
will  be  discussed  and  explained.  More  detailed 
explanations  are  given  in  References  21,  22,  23 
and  24. 

a.  Effects  on  Range 

(1)  Height  of  Target.  Two  methods  may  be 
used  for  determining  die  effect  of  the  height  of 
the  target  with  respect  to  the  gun :  one  by  means 
of  the  trajectory  chart,  and  another  by  a  correction 
to  the  angle  of  departure. 

(a)  For  large  variations  in  height,  the  tra¬ 
jectory  chart  should  be  used.  The  trajectory  that 
passes  through  the  point  determined  by  the  given 
horizontal  range  and  altitude  can  be  found.  The 
total  range  of  this  trajectory  is  the  range  where 
it  crosses  the  axis.  The  difference  between  the 
given  range  and  the  total  range  is  the  effect  of  the 
height  of  target.  For  greater  accuracy,  the  ex¬ 
terior  ballistic  tables  may  be  used  in  lieu  of  the 
charts. 

(b)  For  small  variations  in  height,  an  ap¬ 
proximate  effect  may  be  found  by  adding  the 
angle  of  site, 

«  =  tan-1  (H/R.),  (72) 

to  the  angle  of  departure  obtaining  the  cor¬ 
rected  angle 

f  =  +  (73) 

Here,  H  is  the  height  of  target, 

R, — the  map  range, 

— die  angle  o i  departure  corresponding  to 
the  map  range  in  the  Exterior  Ballistic 
Tables, 

? — the  angle  of  departure  corresponding  to  the 
corrected  range  R. 

Then  the  effect  on  range  doe  to  height  of  target  is 
AhR  =  R.  —  R.  (74) 
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(2)  Elevation.  Since  the  Exterior  Ballistic 
Tables  give  the  range  as  a  function  cf  elevation  far 

constant  values  of  muzzle,  velocity  and  ballistic 
coefficient,  the  effect  cf  an  increase  m  elevation 
may  be  found  from  the  differences.  At  short 
ranges,  the  effect  may  be  compBted  more  ac¬ 
curately  by  using  the  Siacci  tables,  as  explained 
above.  The  vacuum  formula  (Eq,  ^jmay  be  use¬ 
ful  in  estimating  the  effect,  especially  for  lags 
ballistic  coefficients.  It  should  be  noted  that  the 
effect  is  small  in  the  vicinity  cf  the  maximum 
range  and,  cf  course,  is  negative  at  elevations 
above  that  corresponding  to  the  maximum  rargs. 

C3)  Muzzle  Velocity.  Since  the  Exterior  Bal¬ 
listic  Tables  give  the  range  as  a  function  cf  muzzle 
velocity  for  constant  values  cf  elevation  and 
ballistic  coefficient,  the  effect  of  an  increase  in 
muzzle  velocity  may  be  found  from  the  differences. 
At  short  ranges,  the  effect  may  be  computed  by 
using  the  Siacci  tables.  In  vacuum,  the  range  is 
proportional  to  the  square  c£  the  velocity,  as  shown 
in  Equation  6 ;  in  air,  this  is  a  good  approximation 
if  the  ballistic  coefficient  is  large. 

(4)  Ballistic  Coefficient.  Since  the  Exterior 
Ballistic  Tables  give  the  range  as  a  function  cf 
ballistic  coefficient  for  constant  values  cf  muzzle 
velocity  and  elevation;  the  effect  cf  an  increase  in 
log  C  may  be  found  from  the  differences.  At 
short  ranges,  the  effect  may  be  computed  by  using 
the  Siacci  tables.  It  is  customary  to  find  the  effect 
of  a  1  %  increase  in  C,  which  is  equivalent  to  an 
increase  cf  0.00432  in  log  C. 

(5)  Weight  of  Projectile.  An  increase  in 
weight  cf  projectile  decreases  the  muzzle  velocity 
and  increases  the  ballistic  coefficient. 

The  effect  on  muzzle  velocity  may  be  expressed 
ipv0  =nv0(Ap)/p,  (75) 

where  v0  is  the  muzzle  velocity 
p  is  the  projectile  weight, 
n  the  logarithmic  differential  coefficient. 

The  value  of  n  may  be  obtained  frcm  interior 
ballistic  tables.  It  is  approximately  — 0.3  for  a 
rifled  gun  with  multiperforated  propellant  grains, 
—  0.4  for  a  rifled  gun  with  single-perforated 
grains,  —  0.47  for  a  smooth  bore  mortar  with 
flake  propellant,  and  — 0.65  for  a  recoilless  rifle 
with  multiperforated  grains. 

Since  the  ballistic  coefficient  is  proportional  to 
the  projectile  weight  by  dafiniticn,  Equation  (22), 


the  effect  on  C  nay  be  expressed 

A.C=C(Ap)/p.  (76) 

If  A»R  is  the  effect  on  range  cf  1  unit  increase  in 
nuzzle  velocity,  AcR  the  effect  on  range  cf  1% 
increase  m  ballistic  coefficient,  the  effect  on  range 

cf  an  increase  m  projectile  weight  is 

A,R  =  AtR(A,v)  -f  AcR (100  A,C)/C.  (77) 

(6)  Air  Density.  The  retardation  due  to  air 
resistance  is  proportional  to  the  air  density  and 
inversely  proportional  to  the  ballistic  coefficient, 
as  indicated  in  Eq.  (25).  Therefore,  to  a  first 
approximation,  the  effect  on  range  cf  an  increase 
in  air  density  is  equal  to  the  effect  (f  the  same 
percentage  decrease  in  C. 

(7)  Air  Temperature ,  An  increase  in  the 

temperature  cf  the  atmosphere  affects  the  range  in 
two  ways :  it  increases  the  density'  cf  the  air  and 
the  velocity  cf  sound.  The  temperature  cf  the 
air-may  affect  the  temperature  cf  ’the propellant, 
which  affects  the  velocity,  but  that  is  an 

interior  ballistic  problem. 

For  a  given  pressure  and  humidity,  the  density 
is  proportional  to  the  absolute  temperature.  The 
effectof  density  cm  range  has  been  treated  above. 

For  a  given  humidity? the  velocity  cf  sound  is 
proportional  to  the  square  root  cf  the  absolute 
tenperature! .  For  a  given  projectile  velocity,  the 
Mach  number  is  inversely  proportional  to  the 
acoustic  velocity.  A  change  in  M  produces  a 
change  in  the  drag  function.  A  new'  trajectory 
could  be  computed  with  the  varied  function;  how¬ 
ever,  it  is  simpler  to  use  an  approximate  formula 
deduced  frema  semi -empirical  formula  derived  by 
Dr.  GronwaJI,'4* 

A,R  =  0.1929  [0.0IR  -0.005vAR 
—  (1  —  0.000,0l5y,)  A<jR}  At,  (78) 

where  R  is  the  range, 

AXR  —the  increase  in  range  due  to  increase  of 

temperature 

v, — the  muzzle  velocity  in  feet  per  second, 

AtR — the  effect  on  R  cf  1  fps  increase  in  v* 

y* — the  maximum  ordinate  in  feet, 

AcR — the  effect  on  R  of  1%  increase  in  C, 

At — the  increase  in  temperature  in  degrees 
Fanrenbeit 
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(8)  Wind.  A  rear  wind  has  three  effects:  it 
moves  the  coordinate  system  along  with  it  rela¬ 
tive  to  the  earth;  it  increases  the  angle  of  de¬ 
parture  relative  to  the  air;  and  it  decreases  the 
muzzle  velocity  relative  to  die  air.  The  total 
effect  on  range  may  be  expressed  as  follows : 

.  _  nr  Tt  .  sinf  AfR  ArRl 

^=w.  LT+-V-5i— “'id- 

(79) 

where  W*  is  the  range  component  of  die  wind  (a 
rear  wind  is  positive;  a  head  wind  is  negative), 

T — the  time  cf  flight, 
f — the  angle  of  departure, 
v„ — the  muzzle  velocity, 

AfR — the  effect  on  range  due  to  a  change  of 
Af  radians  in  f , 

AtR — the  effect  on  range  due  to  a  change  Av,  in 
v„. 

h.  Effects  on  Deflection 

(1)  Wind.  A  crosswind  has  two  effects:  it 
moves  the  coordinate  system  with  it  and  it  changes 

the  direction  of  motion  relative  to  the  air.  The 
total  effect  on  linear  deflection  may  be  expressed 

A.Z  =  W.  (T  -  R/v.  cos  f ).  (80) 

where  W,  is  the  crosswind. 

The  angular  deflection  due  to  wind,  in  radians, 
is 


D„s=  W,  (T/R—  1/v.cosf).  (81) 


(2)  Drift.  After  the  yaw  due  to  initial  dis¬ 
turbances  has  damped  cut,  a  spuming  shell  has 
a  yaw  of  repose  due  to  the  curvature  cf  the  tra¬ 
jectory.  This  yaw  is  mainly  to  the  right  cf  the 
trajectory  if  the  shell  has  right-handed  spin;  to 
the  left  if  it  has  left-handed  spin.  It-may  be  de¬ 
termined  at  short  ranges  by  firing  alternately 
from  barrels  with  right  and  left  hand  twists  of 
rifling,  or  at  long  ranges  by  correcting  the  ob¬ 
served  deflection  for  the  effects  cf  wind,  cant  cf 
trunnions,  and  rotation  cf  earth. 


The  differential  equation  cf  motion  for  the  drift 

z  is" 


_  2*gAv«LN  .  Gpz 
z  mdnu*MN. X  frC~’ 


(82) 


where  g  is  the  gravitational  acceleration, 

A — the  axial  moment  of  inertia, 
v. — the  muzzle  velocity, 

L — the  cross  wind  force  (lift), 

N — the  gain, 

N, — the  initial  spin, 
m — the  mass  cf  the  projectile, 
d — the  caliber, 

u — the  velocity  cf  the  projectile  relative  to  the 


Mite  overturning  moment  about  the  center  cf 

gravity, 

x — the  horizontal  range, 

G — the  drag  function, 

p — the  air  density, 

f, — the  standard  air  density  (0.07513  lb/ft®), 

C — the  ballistic  coefficient. 

Dots  denote  time  derivatives,  and  the  initial  con¬ 
ditions  are 

z.  =  0,  z.  =0. 

c.  Ballistic  Density,  Temperature,  and  Wind 

In  treating  the  effects  cf  variations  in  the 
atmosphere,  it  has  been  assumed  that  the  percent¬ 
age  variations  from  standard  density  and  tempera¬ 
ture  and  the  wind  are  constant  all  alcng  the  tra- 
jectniY.  Actually,  tb^  vary  with  altitude,  so 
that  weighted  mean  values  must  be  used.  The 
weighting  factors  are  determined  by  the  propor- 
tion  cf  the  total  effect  that  a  variation  occurring 
in  each  altitude  zone  would  have.  Instead  cf 
computing  a  set  of  wed^Ttdng  factors  for  each  tra¬ 
jectory,  a  few  representative  curves  are  used.  The 
sum  cf  the  products  cf  the  zonal  variations  by 
their  weighting  factors  is  called  the  ballistic 
density,  temperature,  or  wind.  The  results  are 
then  treated  as  constants  in  computing  the  effects 
on  range  and  deflection. 
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GLOSSARY 


4  b  of  departure,  (elevation)  As  used  in  this 

pamphlet  these  terms  mean  quadrant  angle  of 
departure  which  is  defined  as  the  angle  between 
the  horizontal  and  the  line  along  which  the  pro¬ 
jectile  leaves  the  projector.  One  cf  the  two 
initial  conditions  cf  a  trajectory,  the  other  con¬ 
dition  being  initial  velocity. 

For  a  rocket-fired  projectile  in  which  the  pro¬ 
pellant  burning  continues  after  the  rocket  leaves 
the  launcher,  it  is  necessary  to  make  assump¬ 
tions  relative  to  the  point  at  which  the  rocket 
becomes  a  missile  in  free  flight  and  the  velocity 
attained  at  that  point. 

ingle  nf  indinatioii.  The  angk  between  the  hori¬ 
zontal  and  the  tangent  to  the  trajectory  at  any 
point.  At  the  origin  cf  the  trajectory  this  is  the 
same  as  the  angle  of  departure. 

tpical  (semi-apical)  angle.  In  general  the  angle 
formed  at  the  apex  or  tip  of  anything.  As  ap¬ 
plied  to  projectiles,  the  angle  between  the  tan¬ 
gents  to  the  curve  outlining  the  contour  cf  the 
projectile  at  its  tip,  or  for  semi-apical  angle, 
the  angle  between  the  axis  and  one  of  the  tan¬ 
gents.  For  a  projectile  having  a  conical  tip, 
the  cone  apex  angle. 

ballistics.  Branch  of  applied  mechanics  which  deals 
with  the  motion  and  behavior  characteristics  of 
missiles,  that  is,  projectiles,  bombs,  rockets, 
guided  missiles,  etc.,  and  of  accompanying  phe¬ 
nomena.  It  can  be  conveniently  divided  into 
three  branches  : — interior  ballistics,  which  deals 
with  the  motion  of  the  projectile  in  the  bore  of 
the  weapon;  exterior  ballistics,  which  deals  with 
the  motion  cf  the  projectile  while  in  flight;  and 
terminal  ballistics,  which  is  concerned  with  the 
effect  and  action  cf  the  projectile  when  it  im¬ 
pacts  or  bursts. 

crosswind  force.  The  component  of  the  force  of 
air  resistance  acting  in  a  direction  perpendicular 
to  the  direction  of  motion,  and  in  the  plane  de¬ 
termined  by  the  tangent  to  the  trajectory  fonned 
by  the  center  cf  gravity  of  a  missile  and  the  axis 
cf  the  missile.  (This  plane  is  called  the  plane 
cf  yaw.) 

drag.  The  component  of  the  force  cf  air  resistance 
acting  in  a  direction  opposite. to  the  motion  cf 
the  center  of  gravity  cf  a  missile. 


elevation.  Sec  angle  of  departure. 

exterior  ballistics.  See  ballistics. 

horizontal  range,  The  horizontal  component  erf 
the  distance  from  the  beginning  of  tie  trajectory 
(point  cf  origin)  to  any  other  point  an  the  tra¬ 
jectory. 

initial  velocity,  (muzzle  velocity)  The  projectile 
velocity  at  the  moment  that  the  projectile  ceases 
to  be  acted  upon  by  propelling  forces. 

For  a  gun -tired  projectile  the  initial  velocity, 
expressed  in  feet  or  meters,  per  second,  is  called 
muscle  velocity  or  initial  velocity.  It  is  obtained 
by  measuring  the  velocity  over  a  distance  for¬ 
ward  d.  the  gun,  and  correcting  back  to  the 
muzzle  for  the  retardation  in  flight. 

For  a  rocket-fired  projectile,  the  velocity  at¬ 
tained  at  the  moment  when  rocket  burnout  oc¬ 
curs  is  the  initial  velocity.  A  slightly  fictitious 
value  is  used,  referred  to  an  assumed  point  as 
origin  at'  the  trajectory. 

The  initial  velocity  cf  a  bomb  dropped  from 
an  airplane  is  the  speed  cf  the  airplane. 

Magaoj  force.  The  component  cf  the  force  cf  air 
resistance  acting  in  a  direction  perpendicular  to 
the  plane  cf  yaw  (see  crosswind  force).  It  is 
caused  by  the  actim  between  the  boundary  layer 
cf  air  rotating  with  the  projectile  and  the  air 
stream.  This  force  is  small,  compared  with 
drag  and  crosswind  force. 

muzzle  velocity.  The  velocity  with  which  the 
projectile  leaves  the  muzzle  cf  the  gun.  See 
initial  velocity. 

slant  range.  The  distance,  measured  along  a 
straight  line,  from  the  beginning  cf  the  trajectory 
(point  cf  origin)  to  a  point  on  the  trajectory. 

stability  factor-  A  factor  which  indicates  the  rela¬ 
tive  stability  of  a  projectile  under  given  condi¬ 
tions.  It  is  dependent  upon  the  force  manents 
acting  to  align  the  projectile  axis  with  the  tan¬ 
gent  to  the  trajectory,  and  the  overturning  force 
moments .  For  the  projectile  to  be  stable,  the 
stability  factor  must  be  greater  than  unity. 

standard  atmosphere.  An  average  or  representa¬ 
tive  structure  of  the  air  used  as  a  point  cf  de¬ 
parture  hi  computing  firing  tables  for  projec- 
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tiles.  The  standard  atmosphere  adopted  in  1956 
for  the  Armed  Services  and  now  known  as  the 
U.  S.  Standard  Atmosphere,  is  that  used  by  the 
International  Civil  Aviation  Organization 
(ICAO) .  ..This  Standard  Atmosphere  assumes 
a  ground  pressure  cf  760  m/m  cf  mercury,  a 
ground  temperature  cf  IS 6  C,  air  density  at 
ground  level  of  0.076,475  Ib/ft4,  and  a  tsrpera- 
ture  gradient  to  the  beginning  cf  the  strato¬ 
sphere  expressed  by  the  fomnila 

absolute  temperature  T(4K)  =  288.16-6.5H 

where  H  is  height  above  sea  level  in  kilometers. 

In  the  stratosphere  (above  11  kilometers)  the 
temperature  is  assumed  constant  at  216,66®  K, 

For  firing  table  purposes  wind  velocities  cf  zero 
are  assumed. 


supercitv.  tion.  An  added  positive  angle  in  anti¬ 
aircraft  gunnery'  that  compensates  for  the  fall  cf 
the  projectile  during  the  time  cf  flight  due  to  the 
pull  cf  gravity;  the  angle  the  gun  or  projector 
nusfc  be  elevated  above  the  gun-target  line. 

trajectory..  The  path  followed  in  space  by  the 
center  of  gravity  of  a  missile  while  in  free  flight. 

yaw.  The  angle  between  the  direction  of  motion 
cf  a  projectile  and  the  axis  cf  the  projectile,  re¬ 
ferred  to  either  as  yaw,  or  more  completely, 
angle  of  yaw,  The  angle  ot  yaw  increases  with 
time  cf  flight  in  an  unstable  projectile  (see  sta¬ 
bility  factor)  and  decreases  to  a  constant  value, 
called  the  yaw  of  repose,  or  the  repose  angle  of 
yaw,  in  a  stable  projectile. 
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ALL  DIMENSIONS  IN  CALIBERS 

Figure  L  Projectile  Shapes 
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Figure  2.  Drag  Coefficient  vs  Mach  Number 


Figure  3.  Time  of  Flight  for  Projectile  Type  1 
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HORIZONH&)  RpJ-oE  (  000 YDS) 

Figure  6  (3  of  12).  Trajectories  for  Proiectile  Tvoe  2.  45*  elevation.  2800  fn«  miiwla  irttl  f«» 
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Figure  6  (8  of  12).  Trajectories  for  Projectile  Type  2,  60'  elevation,  2000  fpe  muzzle  velocity 
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HORIZONTAL  ( GCO  §DS) 

Figure  7  (4  of  12).  Trajectories  for  Projectile  Type  8,  45°  elevation,  3600  fps  muzzle  velocity 
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7  (6  of  12).  Trajectories  Sor 


HORIZONTAL  RANGE  (IQOO  Y08.) 

Figure  7  {T  of  12).  Trajectories  for  Projectile  Type  8,  60'  elevation,  1200  fpr  imussle  velocity 
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Figure  7  (11  of  12).  Trajectories  for  Projectile  Type  8,  60°  elevation,  4400  fpe  muzzle  velocity 
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